Coherent collisional spin dynamics in optical lattices 
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We report on the observation of coherent, purely collisionally driven spin dynamics of neutral atoms in an op- 
tical lattice. For high lattice depths, atom pairs confined to the same lattice site show weakly damped Rabi-type 
oscillations between two-particle Zeeman states of equal magnetization, induced by spin changing collisions. 
This paves the way towards the efficient creation of robust entangled atom pairs in an optical lattice. Moreover, 
measurement of the oscillation frequency allows for precise determination of the spin-changing collisional cou- 
pling strengths, which are directly related to fundamental scattering lengths describing interatomic collisions at 
ultracold temperatures. 

PACS numbers: 03.75.Lm, 03.75.Gg, 03.75.Mn, 34.50.-s 



The creation and manipulation ofspinor Bose-Einstein con- 
densates (BEC) in optical traps 111 yl |4|, |5l |y] has opened a 
wide field of fascinating phenomena originating from the spin 
degree of freedom. Spinor systems allow to explore quantum 
magnetism [7] in a range of parameters not easily accessible 
in solid state physics. They have been proposed e. g. for the 
investigation of spin waves iH 0] or for the study of strongly 
correlated ground states in optical lattices (see lll(illllll2Lll3ll 
and references therein). More importantly, new aspects have 
emerged, such as the generation of entangled states 
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or dynamical s pin mixing induced by interatomic collisions 
Q El OH EllIS HH . The latter phenomenon has been in- 
vestigated in several recent experiments with Bose-Einstein 
condensates in optical dipole traps |3l|4l|5|]. The fundamen- 
tal mechanism responsible for many of these spin phenom- 
ena is a coherent collisional process in which the spin of each 
colliding particle is changed while the total magnetization is 
preserved. 

In this Letter, we investigate this microscopic collision pro- 
cess in an ensemble of isolated atom pairs localized to lat- 
tice sites of a deep optical lattice. We observe coherent os- 
cillations between two-particle Zeeman states, coupled by the 
spin-changing interaction. We show that for a broad range 
of parameters this dynamics can be described by a Rabi-type 
model. Our system allows for a precise measurement of the 
coupling strength for spin changing collisions. Finally, we ob- 
serve an increasing damping of the oscillations with decreas- 
ing lattice depth, which we attribute to the onset of tunneling 
in the system. 

Let us consider a pair of 87 Rb atoms, localized in the vi- 
brational ground state of a deep trapping potential, as shown 
in Fig. [TJ. The trapping frequency oj is assumed to be 
much larger than the typical interaction energy, so that ex- 
citations into higher vibrational levels are suppressed. The 
atom pair can then be solely described by a spin wavefunc- 
tion \fi,muf2,m2), where fi, mi are the total angular mo- 
mentum and its projection onto the z-axis, and i = 1,2 la- 
bels the first and second atom, respectively. In the follow- 
ing, we assume that both atoms are in the upper hyperfine 
ground state with f\ = fi = 2, and abbreviate the non- 
symmetrized two particle states as \m\,rri2). In the absence 



of an external radio-frequency (rf) field, interatomic colli- 
sions drive the spin evolution of this system. For two inter- 
acting alkali atoms, the projection of the total angular mo- 
mentum on the quantization axis is conserved, even in a fi- 
nite magnetic field III |2l yl LUJ ll8l LL3] • The interaction 
thus couples an initial state \<pi) = |mi,m2> to a final state 
= |m3, m4>, provided the total magnetization is conserved, 
i. e. m\ + iri2 - /??3 + m\. Furthermore, s-wave collisions 
between spin f = 2 bosons are characterized by three scat- 
tering lengths of for the collision channels with total spin F 
(F = 0, 2, 4) HI EH EH]. The matrix element Q if of the in- 
teraction hamiltonian between states \<pi) and \<pf) ("coupling 
strength") is proportional to (4nhAa)/M J \(po\ 4 d 3 r, where 0o 
is the (spin-independent) spatial wave function of the ground 
state in the potential well, M is the atomic mass, and Aa is a 
weighted difference of the scattering lengths a? that depends 
on the specific values of the magnetic quantum numbers. For 
example, in the case \<pi) = |0,0) and \<pf) = | + 1,-1), one 
finds Aa = (-lao - 5<22 + 12aA/35. For 87 Rb the different ap 
are almost equal (see e. g. [20, 22]), so that the differences Aa 
are typically only a few percent of the bare scattering lengths. 

Due to the localization in the vibrational ground state and 
due to the constraint of a conserved magnetization, a given ini- 
tial state can couple only to a few final two-particle states par- 
ticipating in the spin evolution (see Fig. [I]}). In the case where 
only one final state is available (relevant for most experiments 
described below), the dynamics reduces to a Rabi-like model. 
We stress, however, that this system is quite different from the 
"usual" Rabi model coupling two single -particle states under 
the influence of an external field. Here we investigate the co- 
herent coupling between two two-particle states induced by 
direct atomic interactions. The Rabi model is parameterized 
by a coupling strength discussed above, and a detuning 
6if = So + S(B 2 ) between the initial and final states, where B 
is the value of the static external magnetic field. As the mag- 
netization is conserved, the initial and final states experience 
the same first order Zeeman shift which has no influence on 
the spin dynamics. However, the second order Zeeman shifts 
are different and introduce the B 2 -dependent detuning. The 
constant detuning Sq originates from the difference in interac- 
tion energies in the initial and final states. The atom pair thus 
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is expected to oscillate between initial and final state at the 
effective Rabi frequency Q^, = yj&ff + 
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FIG. 1 : (a) Two atoms localized in the vibrational ground state of a 
common lattice well can change their spin orientation while preserv- 
ing total magnetization. The atoms remain in the lowest vibrational 
state at all times, (b) The process can be described as a coherent cou- 
pling between two-particle Zeeman states |0 f ), \(pf) and \(f>f>), where 
the coupling strength depends on the choice of initial and final 
state, and the detunings 8 t f, 6ff> can be varied by the second order 
Zeeman shift. 

In our experiment we investigate an array of localized atom 
pairs. We prepare a Mott-insulator of around 2 x 10 5 87 Rb 
atoms in|/= 1, rrif = -l)ina combined optical and magnetic 
trap similar to our previous work 12311 . In this work, unless 
stated otherwise, we use a lattice depth of 40 E r . Here E r = 
h 2 /2MA 2 is the single photon recoil energy, and A « 840 nm 
the lattice laser wavelength. The overall harmonic confine- 
ment of the system leads to the creation of Mott-shells with 
different filling factors (see [23] and references therein). For 
our parameters, we calculate that approximately half of the 
atoms are in the central region with two atoms per site. The 
remaining atoms are distributed in a surrounding shell with 
one atom per site and do not undergo atomic collisions. The 
number of sites with more than two atoms is negligible. 

To study spin dynamics, we subsequently switch the mag- 
netic trap off and load the sample into a pure optical lattice 
lEill . In order to preserve spin polarization of the atoms, a 
homogeneous magnetic field of approximately 1.2 G is main- 
tained B26I1 . The spin dynamics is initialized by transferring 
the sample into either \f - 2, rrif - 0) or \f = 2, rrif = -1), 
and the magnetic field is subsequently ramped to a final value 
between 0.2 G and 2 G. After time evolution for a variable 
time t the optical trap is switched off. In order to spatially 
separate the different magnetic substates, a magnetic gradi- 
ent field is switched on during the first 3 ms of time-of-flight 
(TOF) GLQHEII1. The population^, of each Zeeman level 
rrij is then detected after 7 ms TOF with absorption imaging. 

We first consider the case where we start with both atoms 
in \<pi) = |0, 0). This state couples to \cpf) = | + 1,-1) and 
10^/) = | + 2,-2). The coupling strength for <-> \<j>f,) 
is calculated to be two orders of magnitude smaller than for 



10/), and can be neglected. However, a two-step cou- 
pling channel \<f>i) <-> \<pf) <-> |0y) is also possible, with com- 
parable coupling strengths for each step. Although present 
at low magnetic field |3l|5|], this two-step process is increas- 
ingly suppressed as the magnetic field is increased due to the 
large detuning. For B > 0.6 G, the system mostly oscillates 
between |0, 0) and | + 1,-1). This is shown in Fig. [2| for the 
case of B - 0.8 G. The relative populations in Fig. |2| have 
been calculated as No/N tot for |0, 0) and (N+i + N-i)/N tot for 
1 + 1, -1), where N tot is the independently counted total atom 
number B24I1 . In order to describe the oscillations we use a 
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FIG. 2: Spin dynamics of atom pairs localized in an optical lattice 
at a magnetic field of B = 0.8 G. The atoms are initially prepared 
in |0, 0) and can evolve into | + 1, -1). Shown are the populations 
in m; = (O) and mj - +1 (•) together with a fit to a damped 
sineyielding an oscillation frequency of -2nx 278(3) Hz. 



Rabi-like model, where the transition probability to the final 
state is Q 



Q' 2 f 2 V 

if 



cos(Q^) e^V), 
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with jif being a phenomenological damping rate. The mea- 
sured population in rrif = +1 is simply given by n Pf, where 
n is the fraction of atoms localized in doubly occupied lattice 
sites. A fit to the data using Eq. [fallows to extract the oscilla- 
tion frequency and the damping rate. As explained above, for 
magnetic fields B < 0.6 G, the process | + 1, -1) <-> | + 2, -2) 
starts to play a role. In this magnetic field range we observe 
deviations from the pure sinusoidal evolution, as expected for 
a three level system. Due to limited space, the three level 
model will be presented elsewhere 12511 . 

We have also investigated the case where the sample is pre- 
pared in the initial state |0) = | - 1,-1). The situation is sim- 
pler than in the previous case, because no third level exists 
to which the final state |0, -2) could couple to while preserv- 
ing magnetization. In this case the description as a two level 
system is appropriate at any magnetic field. 

For these two initial two-particle states |0, 0) and | - 1, -1), 
spin oscillations have been observed for various magnetic 
fields, corresponding to different detunings. The measured 
oscillation frequency, plotted in Fig. are well fitted by the 
expected behaviour of the effective Rabi frequency Q!(B) with 
varying detuning (solid lines in Fig.O. As the magnetic field 
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FIG. 3: Oscillation frequency of spin dynamics vs. magnetic 
field (corresponding detuning in upper axis) for the case |0, 0) <-> 
| + 1, -1> (O) and | - 1, -1> <-> |0, -2> (•). The diamonds are re- 
sults of a Fourier transform for the coupled three-level case. The 
solid lines are fits to the expected behaviour of the effective Rabi fre- 
quency. The upper curve has been offset by 300 Hz for clarity. The 
error bars are typically on the order of a few percent. As comparison 
the calculated curves (based on 1 22]) are shown as dashed lines. The 
shaded regions mark the range of frequencies consistent with an er- 
ror of 0.2 a B in the scattering lengths and an uncertainty of 10% (as 
upper bound) in the lattice depth, where both errors have been added. 

dependence S(B 2 ) is known, this fit allows to extract the bare 
coupling strength and zero-field detuning So, related to 
known combinations of the characteristic scattering lengths 
dp in the way explained above. The parameters for the differ- 
ent processes determined by the fits are summarized in Table 
HI These measurements represent a high precision test of spin- 
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243(2) Hz 
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247(3) Hz 
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TABLE I: Summary of measured coupling strength Q if , zero-field 
detuning S and effective Rabi frequency Q! if = ^Jd^ + Q? f at 40 E r 
lattice depth. 

dependent scattering properties and can be used in particular 
to test potentials for ultracold collisions for two Rb atoms. 
Knowledge of these potentials allows to compute the scatter- 
ing lengths a F , from which theoretical values for and So 
can be calculated for our trapping conditions and can be di- 
rectly compared to the experiment. We find a general agree- 
ment in Fig. [3] between the measured oscillation frequencies 
and the calculated ones, based on the prediction of 12211 for 
the ap. However, a discrepancy on the order of 20% remains. 
More details on this comparison will be given in a forthcom- 
ing publication 1 25] . 

The visibility of the oscillations shown in Fig.|2k is smaller 
than the one expected from the Rabi model alone for this de- 
tuning. This is due to the fact that only those atoms with two 
atoms per site contribute to the spin dynamics, whereas the 
relative populations are normalized with respect to the total 



atom number. Hence the visibility is artificially decreased, 
which we accounted for by the factor n introduced above. 
With the extracted values of Q/y and Sif we fit n to the am- 
plitude of the spin oscillation (see solid line in Fig.|4k). The 
mean measured fraction resulting from the fit is n = 0.42(3) 
close to the value n « 50% calculated for the ground state. 

Further evidence comes from a separate experiment, where 
we show that the contrast can be enhanced by selectively dis- 
carding atoms from the measurement which are in sites with 
unity filling. This is accomplished by first evolving the system 
for one half period of spin-oscillation, at which the occupa- 
tion of the | + 1, -l)-state is maximal in sites with two atoms. 
Atoms in the ntf = state, predominantly in singly occupied 
sites, are then transferred into the / = 1 hyperfine manifold by 
the application of a 6 /is microwave 7r-pulse. As the probe light 
of the absorption imaging is resonant with the / = 2 — > /' = 3 
optical transition only, those atoms remain undetected in the 
subsequent spin-evolution. The resulting high contrast oscil- 
lations are shown in Fig.0t for the case of B = 0.6 G. We mea- 
sure an increase of the initial amplitude from 0.31(3) without 
the microwave to 0.56(4). The difference to the theoretically 
expected value of 0.8 can be explained by a small fraction of 
atoms in rrif = ±2 (approximately 17%) 12811 . 
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FIG. 4: (a) Oscillation amplitude vs. magnetic field for the processes 
|00) <-> | + 1, -1) (O) and | - 1, -1) <-> |0, -2) (•). The lines are fits 
to the expected behavior of the Rabi amplitude, where only the frac- 
tion of atoms n in sites with a filling of two remains as free parameter. 
The upper curve has been offset by 0.2. (b) High contrast spin oscil- 
lations between |0, 0) (O) and | + 1, -1) (•) at 0.6 G and 40 E r . Here, 
atoms in sites with unity filling are selectively discarded from the 
measurement (see text). 

The picture of isolated atom pairs coherently oscillating 
between two-particle states does not allow to understand the 
damping of these oscillations, clearly seen in Fig. |5| To in- 
vestigate the damping mechanism, we record the decay rate 
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as a function of lattice depth at B = 0.6 G (see Fig.Q. Here 
the depth has been changed between 20 E r and 54 E r in both 
horizontal lattice axes. For high lattice depths the damping 
rate levels at a finite value, whereas it increases for decreas- 
ing lattice depths. For depths smaller than 20 E r the coher- 
ent oscillations could not be recorded any more. We interpret 
the increase as a signature of the onset of tunneling in the 
lattice. A fit to the damping rate of the form aJ/h + yo re- 
turns a proportionality constant a - 9(1) and a constant offset 
yo = 37(3) s" 1 (see solid line in Fig.|5]for the fit and dashed 
line for the offset). This offset can be explained by the atom 
loss rate in the lattice, as shown in the inset of Fig. For 
this experiment, the sample was prepared in the |0, 0) state 
and spin dynamics was suppressed by a magnetic field around 
10 G. The initial loss rate y\ = 35(5) s" 1 , attributed to two- 
body inelastic processes with a measured two-body loss coef- 
ficient K 2 = (8.8 ± 1.5) x 10" 14 cm 3 /s, agrees with the offset in 
the damping within the quoted uncertainties. We have also de- 
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FIG. 5: Damping rate of spin oscillations vs. lattice depth for |0, 0) <-> 
| + 1,-1>. The solid line is a fit to aJ/h + y with y = 37(3) s _1 
(illustrated as dashed line). The inset shows atom loss vs. hold time 
at a 10 G magnetic field, where spin dynamics is suppressed. The fast 
loss rate of y\ « 35(5) s _1 coincides with the offset yo of the damping 
rate even at high lattice depths. 

tected spin dynamics in the / = 1 hyperfine manifold between 
the two-particle states |0, 0) <-> | + 1 , - 1 ) and find a coupling 
strength roughly one order of magnitude smaller than in / = 2 
I25L in agreement with theoretical predictions. 

In summary, we have observed coherent spin dynamics be- 
tween two-particle states in the upper hyperfine ground state 
of 87 Rb due to spin changing collisions. The observation 
of high contrast Rabi-type oscillations make this system a 
promising starting point for quantum information purposes. 
In this work we have demonstrated a method to create an 
array of entangled atom pairs similar to the long-lived Bell 
pairs produced in ion traps 0> which constitutes a first step 
towards the creation of pair-correlated atomic beams as pro- 
posed in |lJ,Ll5|]. Another intriguing question is the evolution 
of quantum correlations upon melting the Mott-insulator. A 
possible outcome would be a non-local condensate of Bell- 
like pairs delocalized over the entire cloud, which could be 



detected through counting statistics lioll . 
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